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Abstract 

A new version of carlomat that allows to generate automatieally the Monte Carlo programs ded- 
ieated to the deseription of the proeesses e^e~ hadrons at low eenter-of-mass energies is pre¬ 
sented. The program has been substantially modified in order to ineorporate the photon-veetor 
meson mixing terms and to make possible eomputation of the helieity amplitudes involving the 
Feynman interaetion vertiees of new tensor structures, like those predicted by the Resonance Chiral 
Theory or Hidden Local Symmetry model, and the effective Lagrangian of the electromagnetic in¬ 
teraction of the nucleons. Moreover, a number of new options have been introduced in the program 
in order to enable a better control over the effective models implemented. In particular, they offer a 
possibility to determine the dominant production mechanisms of the final state chosen by the user. 

*E-mail: karol.kolodziej @us.edu.pl 



PROGRAM SUMMARY 


Program title: cajcloma.t, version 3.0 
Catalogue identifier: 

Program summary URL: 

Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland 

Licensing provisions: Standard CPC lieenee 

No. of lines in distributed program, including test data, etc.: 

No. of bytes in distributed program, including test data, etc.: 

Distribution format: tar.gz 
Programming language: Fortran 90/95 
Computer: All 
Operating system: Linux 
Classification: 

Nature of problem: 

Predietions for reaetions of low energy -annihilation into final states eontaining pions, kaons, 
light veetor mesons, one or more photons and light fermion pairs within the Standard Model and ef- 
feetive models inspired by the Resonanee Chiral Theory or Hidden Loeal Symmetry model. Deserip- 
tion of the eleetromagnetie produetion of nueleon pairs within the effeetive Lagrangian approaeh. 
Solution method: 

As in former versions, a program for the Monte Carlo (MC) simulation of hadrons at low 

energies is generated in a fully automatie way for a user speeified proeess. However, the user is 
supposed to seleet a number of options and adjust arbitrary parameters in the main part of the MC 
eomputation program in order to obtain possibly the best deseription of experimental data. To this 
end, the user ean also easily supplement her/his own formulae for ^-dependent veetor meson widths 
or running eouplings by appropriately modifying eorresponding subroutines. 

Reasons for new version: 

Proeesses of e^e~ —hadrons in the energy range below the J/yf threshold eannot be deseribed 
in the framework of perturbative quantum ehromodynamies. The sealar eleetrodynamies whieh has 
been implemented in carlomat_2.0 [HI does not provide a satisfaetory deseription either. The most 
promising theoretieal frameworks in this eontext are the Resonanee Chiral Theory or Hidden Loeal 
Symmetry model whieh, among others, involve the photon-veetor meson mixing and a number of 
vertiees of rather eomplieated Lorentz tensor strueture that is not present in the Standard Model or 
sealar QED. Already at low energies, the hadronie final states may eonsist of several partieles, sueh 
as pions, kaons, or nueleons whieh ean be aeeompanied by one or more photons, or light fermion 
pairs sueh as e^e^, or p^pL. The number of Feynman diagrams of sueh multipartiele reaetions 
grows substantially with inereasing numbers of interaetion vertiees and mixing terms of the effee¬ 
tive models. Therefore, it is highly desirable to automatize the ealeulations. At the same time, new 
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program options should provide the user with an easy way of implementing her/his own changes in 
the program in order to better fit the experimental data. 

Summary of revisions: 

The code-generation part of the program has been substantially modified in order to incorporate the 
photon-vector meson mixing and calls to new subroutines for computation of the helicity amplitudes 
of the building blocks and complete Feynman diagrams which contain new interaction vertices and 
mixing terms. The subroutine library of carlomat has been extended to make possible computation 
of the helicity amplitudes involving the Feynman interaction vertices of new Lorentz tensor struc¬ 
tures. Many subroutines have been modified in order to incorporate the ^^-dependent couplings and 
vector meson widths. A number of options have been introduced in order to give a better control of 
the effective model implemented. 

Restrictions: 

As in previous versions of the program the number of particles is limited to 12 which exceeds typical 
numbers of particles of the exclusive low energy -annihilation processes. However, in the pres¬ 
ence of photon-vector meson mixing, the Feynman diagrams proliferate, for example, with currently 
implemented Feynman rules, there are 90672 diagrams of —)■ 3(7t+7t^). Hence, the compi¬ 

lation time of generated code may become very long already for processes with smaller number of 
the final state particles. Many couplings of the effective models are not known with good enough 
precision and must be adjusted in consecutive runs of the program in order to obtain satisfactory 
description of the experimental data. 

Running time: 

Depends on the selected process. Typical running time for the code generation vary from a fraction 
of asecondfor, e.g., —)■ to about 2 minutes for —)■ 3(7t+7t^). It may become 

substantially longer for processes with more particles in the final state. The execution time necessary 
to produce the appended test output files for e^e~ TtfTtrg&gry and e^e^ —)■ 7t+7t^7t+7t^Y was 
13s and 4s, respectively. The code generation for both processes took a fraction of a second time for 
each process. 


1 Introduction 


Hadronic contributions to the vacuum polarization are the major factor that influences precision of 
theoretical predictions for the muon anomaly and plays an important role in the evolution of the 
fine structure constant a(2^) from the Thomson limit to high energy scales. Improving the precision 
of predictions for the muon anomaly becomes vital in the prospect of forthcoming measurements in 
Fermilab that should reduce the experimental error of to 0.14 parts per million, while the better 
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precision of would be important for the precision data analysis from the future high energy 

collider, which would most probably include a giga-Z option. Because of the breakdown of 
predictive power of the perturbative QCD at low momentum transfer, the hadronic contributions to 
the vacuum polarization are determined, with the help of dispersion relations, from the energy de¬ 
pendence of the total cross section of electron-positron annihilation into hadrons, Og+g-^hadrons('^)- 
Below the J/\\f production threshold, Og+g-^hadrons niust be measured and confronted with theoreti¬ 
cal predictions of some effective model for the low energy hadron physics. 

There are two QCD inspired theoretical frameworks which seem to be applicable in this context: the 
Resonance Chiral Theory (RxT) |I21 and the Hidden Local Symmetry (HLS) model O, which were 
proven to be essentially equivalent [j4l|. For example, the HLS model allowed for a quite satisfac¬ 
tory simultaneous description of most hadronic -annihilation channels in the low energy range, 
including (|)-resonance and 10 decay widths, mostly radiative ones, of light mesons and allowed 
to resolve the inconsistency between the e+e”-annihilation to 7t+7t^ and the x'^-decay to 
[|5l , dUl. The hadronic currents based on R^T were implemented in TAUOLA, a x-decay Monte 
Carlo (MC) generator ||71, and used for description of the x lepton decay into two or three pseu¬ 
doscalar mesons that constitute 88% of the x hadronic decay width in Ref. [[8]| and later improved for 
x"^ —)■ decay mode in Ref. flUl which allowed to successfully model the one-dimensional 

distributions measured by the BaBar collaboration. 

The number of Feynman diagrams in the framework of RxT or HLS model grows quite fast with 
the number of particles in the final state of e~^e~ —?■ hadrons. In particular, in the presence of 
one or a few photon-vector meson mixing terms, it can easily reach a hundred thousand already 
for 67t. Obviously, preparation of a reliable MC generator for such a process is rather 

tedious a task, unless the process of code writing is fully automatized. The first step toward the 
automatic generation of the MC programs for description of reactions —)■ hadrons was already 

made in carlomat_2.0 m, in which the Feynman rules of the scalar electrodynamics (sQED) were 
implemented in addition to those of the Standard Model (SM). It allowed to effectively describe 
the electromagnetic (EM) interaction of charged pions which, at low energies, can be considered as 
being point like particles, see, e.g., ifT^ . However, the charged pion form factor that would allow 
to account for the bound state nature of the particle was not implemented in carlomat_2.0. In the 
present paper, a new version, labeled with 3.0, of a program carlomat [fTTll . [fT]| is described, which 
to large extent should meet the requirements of automatic code generation for MC simulation of the 
low energetic e+e -annihilation into hadrons in the framework of the effective models. 

The paper is organized in the following way. In Section 2, the implementation of new Eeynman 
rules in the program is described. New options which have been implemented in the program to 
give the user a better control over the model are described in Section 3. Einally, the instructions for 
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preparation for running and usage of the program are given in Section 4. 


2 New Feynman rules implemented in the program 


In this section, the implementation of the Feynman rules of the HLS model that are relevant for the 
description of hadrons in the low energy range in carlomat_3.0 is described. Most of the 

rules can be derived from the Lagrangian pieces of Appendix C of Ref. |[5l . The Lagrangian of EM 
interaction of spin 1/2 nucleons implemented in the program is described in Subsection 2.3. 


2.1 Photon-vector meson mixing 

The topology generator of carlomat takes into account only triple and quartic vertices, therefore 
the mixing should be added in subroutine checktop, where topologies of diagrams are confronted 
with the implemented Feynman rules. The procedure was described in detail in Ref. [fTTll . For 
the sake of clarity let us remind here, that every topology in carlomat is divided into two parts, 
each being checked against the Feynman rules separately. This is done by consecutive calls to 
subroutine genpart that combines two (three) particles into the third (fourth) leg of a triple (quartic) 
Feynman vertex which is then folded with the adjacent Feynman propagator to form an off-shell 
particle. The latter is represented by an array of spinors, polarization vectors or scalars, whose 
elements are labeled with different combinations of the polarization indices of the particle spinors or 
polarization vectors of which they are formed. At this point, if the particle mixing is present, a new 
subroutine mixpl is called to check whether the propagator of the off-shell particle can be mixed 
with some other propagator or not, if so, a new off-shell particle is formed. The particle with mixing 
is appropriately tagged in order not to be mixed again, because, according to Fig. [H the mixing 
term contains an extra power of the electric charge e and therefore should be considered as the next 
to leading-order correction. This procedure is being repeated until finally two off-shell particles 
are formed, corresponding to both parts of the considered topology. Then, a modified subroutine 
matchkk is called which checks, whether the particles can be matched in the Feynman diagram with 
the Feynman propagator or, if none of them has been mixed yet, with some of the mixing terms of 

Fig.m 

In spite of being conceptually quite simple, the implementation of particle mixing required sub¬ 
stantial changes in the code-generation part of the program. Moreover, new subroutines bbkk and 
bbmd have been written to compute, respectively, the polarization vectors of the off-shell particle and 
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helicity amplitudes of the Feynman diagrams in the case of mixing. 


A^‘(q) V^q) 

'V\AA^ = -efAv{q^)g‘", Mh V = p“,0),i|i,pi,p2. 

Figure 1: The photon-vector meson mixing diagrams implemented in the current version of the program; pi 
and p 2 stand for p(1450) and p(1700), respectively. 


2.2 Interaction vertices 

The triple and quartic interaction vertices of the HLS model that are implemented in carlomat_3.0 
are depicted in Figs. [21-111 where all the particle four momenta are assumed to be incoming to the 
vertex and with = 1, is the totally antisymmetric Levi-Civita tensor. A number of new 
subroutines for computation of the building blocks and complete amplitudes of the Feynman dia¬ 
grams containing vertices of new tensor structure have been written. The implementation of calls 
to the new subroutines required some changes in the code-generation part of the program which 
concerned mainly subroutine genpart. 

The triple vertices of the photon or vector meson interaction with the pseudoscalar meson pair 
PP which have the form similar to the triple vertex of sQED are shown in Fig.|2l The only difference 
is the replacement 

e efApp{q^) and e fvpp{q^), (1) 

where V = p®,a),(|),pi,p 2 and P = Although couplings of pi = p(1450) and p 2 = 

p( 1700) to other particles are hard to define on the basis of existing data lIT^ . the interaction vertices 
p,7r+7t^ and mixing terms y— p/, i = 1,2, have been included in the program just to enable tests of 
their possible influence on some observables, e.g., on the pion form factor, where they play a role. 

Subroutines ppakk, appkk and papkk for the computation of building blocks of the Feynman dia¬ 
grams, and ppamd, appmd and papmd for the computation of the helicity amplitudes in the sQED of 
carlomat_2.0 have been all supplemented with an option ig, which allows to take into account the 
y^-dependent couplings of ([U). Subroutines ppakk and ppamd have been additionally supplied with 
an option iwdth, that gives a possibility to include the ^'-dependent width of a vector meson. The 
use of both options is explained in Section 3. 

Triple interaction vertices of the HES model that have a form different from that of the triple vertices 
of the SM or sQED are depicted in Eig. [3l New subroutines that have been written in order to 
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Figure 2: Triple vertices of the photon or vector meson V = p**, co, (j), interaction with the pseudoscalar 
meson pair PP, P = of the same form as that of the triple vertex of sQED. 


Af^{pi) 


Af^{pi) 




Figure 3: Triple vertices of the pion interaction with photons or vector mesons, where, in the top right corner, 
V = p°, CO, in the bottom right comer P = 11 ^ and T = p*’ or P = 7t^ and T = p^. 

compute the corresponding building blocks and helicity amplitudes are: pvvkk, pvvmd, vvpkk and 
vvpmd. All of them include the running-coupling option ig, and the first one, whose output is an 
array of four vectors corresponding to all possible helicities of the scalar and vector particles they 
are composed of, includes in addition the running-width option iwdth. 

The quartic interaction vertices of the HLS model implemented in the current version of the program 
are shown in Fig. |4l The vertices in the first row have the same tensor form as the quartic vertex of 
the sQED or the quartic vertices of the Nambu-Goldstone boson - gauge boson interaction of the 
SM, which have been implemented already in the first version of carlomat. Hence, the correspond¬ 
ing building blocks and helicity amplitudes can be computed with modified subroutines vvsskk, 





Figure 4: Quartic vertices of the HLS model implemented in the current version of the program. The quartic 
vertex of sQED, implemented already in carlomat_2.0, is not shown. 


vvssmd, vsvskk, vsvsmd, vssvkk, vssvmd, svvskk, svvsmd, ssvvkk, ssvvmd, svsvkk 
and svsvmd, which have been all supplied with the running-coupling option ig. Subroutines vssvkk, 
ssvvkk and svsvkk have been moreover supplemented with the running-width option iwidth. The 
tensor form of the vertices in the second row of Fig. |4] is different. Therefore, the correspond¬ 
ing building blocks and helicity amplitudes are computed with newly written subroutines pppvkk, 
pppvmd, vpppkk and vpppmd. 


2.3 Electromagnetic interaction of nucleons 


The Lagrangian of EM interaction of spin 1/2 nucleons has the following form: 


^ANN = eA/jN (p) 


fFi{Q^) + ^of^^q^F2{Q 


N{p), 


( 2 ) 


where o^'' = ^[^',7''], with y", /r = 0,1,2,3, being the Dirac matrices, q = p — p' is the four mo¬ 
mentum transfer, F\{^) and F 2 {Q^) are the form factors and = —q^. The form of Eq. dH) is 
similar to that of the effective Lagrangian of the Wtb interaction given by Eq. (4) of Ref. [|T||. Due 
to this fact, the implementation of the corresponding Feynman rules for the nucleon-photon interac¬ 
tion was straightforward. To compute the corresponding building blocks and helicity amplitudes the 
following new subroutines have been written: annkk, annmd, nnakk, nnamd, nankk and nanmd 
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and the calls to them have been appropriately implemented in subroutine genpart. The form factors 
^i(2^) and F 2 {Q^) have been adopted from PHOKARA IfTJl with the help of an interface subroutine 
nuclff_phok. In this way, the MC simulations of processes involving the EM production of the 
nucleon pairs have become possible. 


3 New program options 


New options which have been added in the program to give the user a better control over the imple¬ 
mented models for the description of the electron-positron annihilation into hadrons at low energies 
are explained below. 

All subroutines that are used to compute the building blocks or the complete helicity amplitudes of 
the Feynman diagrams of Figs. [21-111 have been supplied with the running-coupling option the name 
of which is formed by adding a prefix i to the name of the corresponding coupling, as the name 
is created in exactly the same way at the stage of code generation. The options are to be specified 
in subroutine couplsm, where they are defined below the assignment instruction for each particular 
coupling. 

icoupl_name=0/l, 2, . . . if the f ixed/running coupling is to be used in the computation, 

where choices 1,2,... corresponding to different running couplings / of Figs. HUH should be 

added by the user as extra else if (ig == . . .) then blocks in subroutine runcoupl. The block 
must contain an assignment for a double complex variable rg in terms of the four momentum transfer 
squared and any other physical parameters that are available in module inprms. The actual form 
of the four momentum transfer q is determined automatically from the four momentum conservation 
in the corresponding interaction vertex at the stage of code generation. Many couplings of the RxT 
or HFS model are not known well enough and therefore must be adjusted in consecutive runs of the 
program in order to obtain satisfactory description of the experimental data. If there are no hints 
as to the form of the running couplings /...(<?^) then it is recommended to set the corresponding 
running-coupling option to 0, which means that the fixed coupling is to be used in the computation. 
The user can also modify any of the fixed couplings by changing the corresponding assignments in 
couplsm, where the couplings are defined in terms of the physical parameters of module inprms. 

The subroutines for computation of the four vectors representing vector mesons have been in addition 
supplied with the running-width option iwdth_name, i.e. igmrh, igmom, igmph, igmrl, igmr2 
for the running width of p°, (0,(|), pi, P 2 , respectively: 
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iwdth_name=0/l ,2,3 if the f ixed/running width of the veetor partiele should be used, 

where ehoiees 1,2,3 refer to different running-width options in subroutine runwidth whieh again 
ean easily be extended by the user. The options are eontrolled from carlocom, the main part of the 
MC eomputation program. 

The main part of the MC eomputation program carlocom contains a few flags: iarho, iaome, 
iaphi, iarhol and iarho2 that allow to switch off and on the photon mixing with p,(0, (|), Pi 
and p 2 vector mesons without a need of running the code-generation program anew, provided that 
the corresponding mixing terms were included in a file vertices.dat when the MC code was 
generated. This gives a possibility to determine the dominant production mechanisms of the final 
state considered by the user. 

In order to give a better control over the mixing contributions to a given process, subroutines bbkk 
and bbmd are equipped with the option: 

iwgt=0/l,2, . . . if the additional complex factor Cl, C 2 , ... is not/is to be included in/ av(^^) 
of Fig. [U i.e. in the amplitude of the Feynman diagrams containing this 
particular particle mixing contribution. 

The actual names for that option in carlocom are: imrho, imome, imphi, imrhl, imrh2 for 
the p®, CO, (|), pi, p 2 meson, respectively. The complex factor cj, j = 1,2,... is given by 

Cj = Wj (3) 

where wj is a positive weight, cpj is an angle in degrees, which should be both specified for each 
possible particle mixing term in the main program for the MC computation carlocom, and fj{q^) 
is a possible four momentum transfer dependence that is defined in subroutine weight!actor. 
Actually only three simple dependences corresponding to iwgt=l,2,3 are currently defined in 
weight!actor, but the user can easily add more options by implementing new else i! (iwgt 
== . . .) then conditions. 

An important new option in the program, which allows to test the EM gauge invariance for processes 
with one or more external photons, is igauge in carlocom.!: 

igauge=l, 2, . . . /else if the gauge invariance is/is not to be tested, 

where 1,2,... is the number of a photon, counting from left to right, whose polarization four vector 
is replaced with its four momentum. 
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To illustrate how this option can be used in practice, consider the following radiative processes: 

e~ —)■ (4) 

-f 7t^7t^7t+7t^y. (5) 

Taking into account the Feynman rules of SM, without the Higgs couplings to electrons and muons, 
sQED, the y— p® mixing of Fig. [Hand the vertices: ytt+Tt^ and of Fig. |2l Tt^yy and 7 r°yp® 

of Fig. |3] and yp® 7 t+ 7 r^ and yjt® 7 t+ 7 r^ of Fig. m processes (jH) and ([5]) receive contribution from, 
respectively, 209 and 774 Feynman diagrams. If, in addition, the vertices Tt^yp"^ of Fig. |3] are 
included then the number of diagrams of processes (jH) and ([5]) grows, respectively, to 231 and 968. 
The cross sections of processes dH) and dU) at 1/5 = 1 GeV, with the following cuts on the angles 
between the photon and a lepton Gy/, the photon and a pion B^i and the photon energy: 

0y/ > 5°, Gy;! > 5°, Ey > 10 McV, ( 6 ) 

are presented in Table 1. The cross sections without (with) contribution from the Tt^yp^ interaction 
vertices of Fig. [3] are printed in the first (second) column for each process. If igauge=l then the 
cross section drops by about 32 orders of magnitude, which means that the EM gauge invariance 
works perfectly well. However, if the vertices Tt^yp^ of Eig. |3] are included then the EM gauge 
invariance is not so perfect any more. Eor process dl]) this is caused by the two Eeynman diagrams 
depicted in Eig. |5l To justify this statement, let us denote the four momenta of particles of process dl]) 
by pi,P 2 ,...,P 7 , from left to right consecutively, and consider the EM gauge invariance test for the 
amplitudes of the diagrams (a) and (b) of Eig. [51 which means in practice that the photon polarization 
four vector is replaced with its four momentum. Neglecting the i factors, which are the same for both 
amplitudes, and skipping polarization indices the amplitudes read: 

I g/'gp 

2 

= ^2 ^j^2 + Pl)^, ( 7 ) 

Mb = g^ei2ve'^^^*^Pl2(x(-^|3) ^2^^^' e56ge^^^^(-P56y)^6‘^47 

2 

= - ^2 ^2 e 12V P12aP 3 [3 £56 ^56 (P 4 + P 7 )^, (8) 

where £j 2 (e^g) is the polarization four vector representing the e+e^y vertex contracted with 

the adjacent photon propagator, = nip — inipYp is the complex p meson mass parameter, pn = 
Pi + P 2 , P 56 = P 5 +P 6 , y = P 56 + P 3 + Pi and r = p 5 (, +P 4 +P'] are four momenta of intermediate 
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virtual photons and mesons, and the eoupling g = = efT^+Ap-{q^) = efn-Ap+{f^) = 

efn+Ap-{f^) has been assumed to have a fixed value. In the seeond row of Eqs. ([7]) and ([ 8 ]), use has 
been made of the faet that, in the EM gauge invarianee test, the sealars 537 and 547 representing the 
and vertex multiplied with the adjaeent pion propagator, take the following form: 


S 31 

547 


(p3+P7-(-P3))^e^(P7) 

{P3+P7)^-ml 

(-P4-(P4+P7))^e*(p7) 

{p4+P7)^-ml 


£{pi)->Pi 

e(P7)->P7 


^2p3_^ 
2p3 ■ pq 

^ ^PA■Pi 
IpA-Pl 


It is elear from the form of Eqs. (|7]) and ([ 8 ]) that amplitudes Ma and Mq, neither vanish separately nor 
eaneel eaeh other, eontrary to the amplitudes of the other 20 Eeynman diagrams of proeess dH) whieh 
also eontain the vertiees Although that degree of gauge invarianee violation should not play 

any role in praetiee, sueh effeets should be treated with great eare, as they may beeome sizable in 
some regions of the photon phase spaee. Therefore, it is reeommended to use the igauge option 
whenever new interaetion vertiees are added to the program. 




Eigure 5: The Feynman diagrams of process ([Hi that spoil the EM gauge invariance. The blobs indicate the 
vertices Ti^yp^. 


4 Preparation for running and program usage 


carlomat_3.0 is distributed as a single tar.gz arehive carlomat_3.0.tgz whieh ean be down¬ 
loaded from: http://kk.us.edu.pl/oarlomat.html. When untared with a eommand 
tar -xzvf carlomat_3.0.tgz 
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igauge 

o(e+e ^ 7t+7t 7t+7t y) 

(5{e^e -> 7t+7t y) 

0 

11.86(5) 11.83(5) 

0.0590(2) 0.0586(2) 

1 

0.124(2)e-30 0.441(l)e-10 

0.636(9)e-33 0.973(l)e-9 


Table 1: The cross sections in pb of processes (011 and (l5]l at ^/s = 1 GeV without (first column) and with 
(second column) contributions from the Tt^yp^ interaction vertices of Fig. [3] The cuts used in the computation 
are given by ([Q. The numbers in parentheses show the MC uncertainty of the last decimal. 


it will create directory carlomat_3.0 with sub directories: code_generation, nic_computation, 
carlolib, test_output and test_outputO. 

Although carlomat_3.0 is dedicated to the description of low energy scattering, interfaces 
to the parton density functions are kept. Therefore, files mstwpdf . f of MSTW lfT4l and Ctq6Pdf . f 
and cteq61.tbl of CTEQ6 ifTSl are also included in the current distribution of the program, but 
grids must be downloaded from the web page of MSTW, see readme file or [|T1| for details. If the 
program will not be run for the hadron scattering processes the user can comment lines contained 
between ckk_had> and ckk_had< in crosskk.f and parfixkk.f, and comment or remove ref¬ 
erences to mstw_interf ace. o, mstwpdf. o, CtqGPdf.o, ctq6f_interf ace. o from makefile 
in mc_computation. 

Preparation for running requires basically the same steps as in carlomat_2.0. They are recollected 
below for user’s convenience. 


• Choose a Fortran 90 compiler in makefile’s of code_generation and mc_computation and 
compile all the routines of carlolib with the same compiler as that chosen in mc_computation; 

• Specify the process and required options in carlomat. f and execute make code from the 
command line in code_generation; 

• Go to mc_computation, choose the center of mass energy and required options in carlocom. f 
and execute make me in the command line. 


Whenever the Fortran compiler is changed, or a compiled program is transferred to another com¬ 
puter with a different processor, all the object and module files should be deleted by executing the 
commands: 
rm *. o 
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rm *.mod 

and the necessary steps of those listed above should be repeated. 

The basic output of the MC run is written to file tot_name, where name is created automatically if 

the assignment for character variable 

prcsnm='auto’ 

in carlomat. f is not changed to arbitrary user’s defined name. The output files for processes dH) 
and ([5]) with the preselected parameters and options should reproduce those delivered in directory 
test_outputO. 

If the differential cross sections/distributions are required then set 
idis=l 

in carlocom. f. The number of distributions to be calculated must be specified in distribs . f and 
their parameters should be defined in calcdis . f. The output will be stored in data files db#_name 
and dl#_name which can be plotted with boxes and lines, respectively, with the use of gnuplot. 
When the run is finished all output files, except for test that may contain information relevant in 
case of unexpected program stop, are moved to directory test_output. 

As in former versions of the program, there is a possibility of generating the unweighted events. It 
is governed by the option imc that is available in carlocom. 

The code generation for processes (01) and Q takes a fraction of a second time. The MC computation 
of the cross sections of Table [3] in 10 iterations, with a maximum of 200000 calls to the integrand 
each, takes 142s and 43s time, respectively, for (jH) and dS]) on processor Intel'’’ Core™ i5-4200M 
CPU @2.50 GHz with a 64 bit Intel Fortran compiler. 
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